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ABSTRACT 


The  wind  stress  calculated  from  wind  velocities  measured  during 
the  Mixed  Layer  Experiment  (MILE)  was  used  as  input  to  a  one-dimensional 
wind-driven  current  model.  These  model  results  are  compared  with  ob¬ 
served  currents  from  the  MILE-1  bouy,  showing  a  qualitative  agreement. 

MILE  was  an  examination  of  the  upper  ocean  carried  out  near  Ocean 
Weather  Station  Papa  during  a  20-day  period,  August  and  September  of 
1977,  which  was  characterized  by  two  major  wind  events. 

The  observed  currents  have  been  analyzed  to  obtain  information  about 
their  behaviour  that  could  be  used  in  the  tuning  of  the  model.  For  a 
simulation  of  the  entire  period  the  results  are  considered  only 
satisfactory. 
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I.  INTRODUCTION 


A.  PURPOSE  OF  THE  STUDY 

The  objective  of  this  research  was  to  investigate  how  well  an  one¬ 
dimensional  model  of  the  wind  driven  current  could  explain  the  observed 
near  surface  current  during  the  Mixed  Layer  Experiment  (MILE).  Measure¬ 
ments  of  the  near-surface  currents  were  taken  during  MILE,  August- 
September  1977,  in  the  vicinity  of  Ocean  Weather  Station  (OWS)  PAPA. 

During  this  period,  the  data  show  the  response  of  the  upper  ocean  to  the 
passage  of  two  wind  events. 

This  study  investigates  the  hypothesis  that  surface  wind  can  account 
for  a  large  part  of  the  energy  and  variability  of  inertial  oscillations 
near  the  ocean  surface  and  that  these  oscillations  are  predominantly 
locally  generated,  as  was  demonstrated  by  Pollard  and  Millard  (1970)  at 
another  site. 

To  accomplish  this  task,  a  one-dimensional  wind-driven  current  model 
was  solved  numerically  using  a  modified  "leapfrog"  scheme.  The  modifi¬ 
cation  involved  employed  current  values,  u  and  v,  which  were  vertically 
averaged  at  each  time  step  to  a  depth  equal  to  that  of  the  observed  mixed 
layer,  simulating  strong  boundary  layer  mixing.  The  mixing  depth  used 
was  a  time-average  of  the  mixed  layer  depth  obtained  from  buoy  temperature 
records  and  verified  with  Plessey  CTD  data  taken  from  the  R/V 
OCEANOGRAPHER. 

The  data  used  to  check  the  response  of  the  model  were  those  obtained 
by  nine  near-surface  current  meters,  at  depths  from  5  to  32m,  of  the 
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MILE-l  flooring.  The  wind  stress  was  computed  from  wind  speed  and  di¬ 
rection  as  recorded  by  the  R/V  OCEANOGRAPHER  during  the  experiment. 

It  was  hoped  that  the  non-wind-driven  current  would  be  sufficiently 
small  and  invariant  so  that  the  hypothesis  could  be  tested. 

In  the  next  section,  the  current  meter  observations  are  analyzed  by 
using  a  Fast  Fourier  Transform.  The  power  spectra  of  currents  from  three 
significant  depths  illustrate  how  the  currents  are  dependent  on  the 
surface  wind  stress. 

In  section  three  the  equations  on  which  the  model  was  based  are  given. 
It  is  shown  how  from  those  equations  a  numerical  solution  is  obtained  to 
solve  the  initial  value  problem  that  we  are  dealing  with  and  how  the 
model  results  were  computed. 

In  the  last  section  the  results  and  conclusions  are  given. 

If  there  is  no  other  indication,  the  abscissa  in  all  figures  gives 
time  in  hours  from  the  start  of  the  experiment,  0500  GI'IT  hours  August 
18th.  There  are  460  hours  of  recorded  currents  in  all. 

3.  IMPORTANCE  OF  THE  PROBLEM 

As  Csanady  [1981]  explains:  "The  problem  of  circulation  is  to  de¬ 
scribe  and  understand  the  pattern  of  the  residual  or  longer-term  water 
particle  displacement.  The  distribution  of  important  water  properties, 
e.g.,  temperature,  salinity  (...),  depends  critically  on  the  pattern  of 
circulation."  The  capability  of  modeling  wind-driven  currents  could 
help  to  explain  one  component  of  that  pattern  of  circulation. 

Modeling  wind-driven  current  is  needed  in  order  to  improve  upper  ocean 

thermal  structure  models,  which  have  great  importance  in  ASW.  Fisheries, 

search  and  rescue  operations,  and  control  of  pollution  in  the  ocean  are 

other  fields  where  wind-driven  currents  must  be  considered. 
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C.  LITERATURE  REVIEW 

Horizontal  inertial  currents  occur  in  thin  layers  with  a  relatively 
reduced  horizontal  extent,  Webster  [1968]  concludes  that  inertial 
oscillations  are  essentially  transient  phenomena.  For  a  fixed  location 
inertial  oscillations  at  different  depths  may  have  diverse  origins. 
According  to  Pollard  and  Millard  [1970],  Gonella  [1971]  and  Kundu  [1976], 
for  the  surface  layer  of  the  ocean,  the  wind  stress  is  the  principal 
mechanism  for  their  generation. 

Inertial  currents  rotate  clockwise  (Northern  Hemisphere)  with  a 
period  T  =  Tr/(li^l  sin  9),  where  is  the  earth's  rotation  vector  and 
9  is  the  latitude  of  the  observation  point.  The  transmission  of  the 
movement  from  the  surface  to  deeper  layers  may  normally  be  carried  out 
either  by  turbulence  in  an  unbounded  homogeneous  medium,  or  by  boundary 
effects  in  a  strongly  stratified  medium. 

To  test  the  dependence  of  inertial  oscillations  as  a  function  of 
wind  stress  the  theory  of  Ekman  [1905]  was  applied  in  an  impulsive 
system.  In  this  theory,  the  wind  provides  a  stress  on  the  sea  surface, 
and  the  motion  is  viscously  transmitted  to  lower  layers. 

Figures  1  and  2  show  four  typical  examples  of  the  measured  vertical 
structures.  In  figure  1  both  temperature  and  velocity  profiles  are  con¬ 
tinuously  stratified:  the  negative  gradients  of  temperature  show  that 
there  is  no  mixing.  However,  the  profiles  in  figure  2  show  a  change  in 
temperature  of  less  than  0.2°  C  in  the  upper  28m,  indicating  that  the 
water  is  nearly  homogeneous  and  the  current  is  also  nearly  uniform.  The 
more  the  ocean  is  stratified,  the  greater  is  the  tendency  for  agreement 
with  Ekman 's  model  which  was  developed  for  a  homogeneous  ocean.  As 
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Gonella  [1971]  concludes  from  his  work  in  the  Mediterranean  this  agree¬ 
ment  comes  about  "because  the  assumption  of  a  constant  coefficient  of 
eddy  viscosity  is  more  realistic  in  a  continuously  stratified  ocean  than 
in  an  entirely  homogeneous  ocean  where  the  turbulence  may  be  too  intense" 

The  approach  to  the  problem  in  this  study  is  similar  to  the  “slab" 
model  of  Pollard  and  Millard  [1970]  and  Garwood  [1976],  using  a  depth 
integral  of  the  momentum  equation  to  solve  the  problem  of  turbulent 
momentum  transfer  within  the  boundary  layer. 

The  essential  principle  is  that  the  integrated  inertial  force,  in¬ 
cluding  Coriolis  and  local  time  derivative  accelerations,  just  balances 
the  applied  surface  stress  less  a  small  damping  term. 

To  interpret  actual  current  records,  two  modifications  were  intro¬ 
duced  in  the  present  work:  (i)  the  momentum  equation  was  integrated  to 
a  depth  equal  to  that  of  the  mixed  layer  in  each  time  step,  and  (ii)  a 
linear  damping  term  was  applied.  By  the  former  procedure,  the  modeled 
mixed  layer  becomes  unstratified  with  f  =  Ztt/T  as  its  only  natural  fre¬ 
quency.  The  latter  step  has  the  effect  of  changing  that  frequency  to  a 
slightly  smaller  value.  This  may  not  be  a  good  property  when  modeling 
the  ocean,  where  inertio-gravity  waves  have  frequencies  between  f  and 
N,  the  Brunt-Vaisala  frequency.  Since  N  is  usually  greater  than  f, 
the  observed  near-inertial  currents  frequently  have  a  frequency  slightly 
greater  than  f,  [Pollard  and  Millard,  1970]. 

In  the  one-dimensional  model  used  here  surface  winds  are  expected  to 
account  for  a  large  part  of  the  variability  of  the  energy  of  the  near- 
inertial  oscillations  near  the  surface.  It  was  also  expected  that,  as 
in  other  models  [Pollard  and  Millard,  1970;  Kundu,  1976],  the  changes  in 


phase  and  amplitude  of  the  forced  oscillation  could  be  simulated  with 


some  success. 

Another  problem  to  be  considered  is  how  the  energy  at  inertial  fre¬ 
quencies  is  removed  from  the  near-surface  layers.  Pollard  [1970], 

Kroll  [1975],  and  Kase  [1979]  show  that  one  possibility  is  that  energy 
is  dispersed  through  the  seasonal  thermocline  by  transferring  into  in¬ 
ternal  waves  in  the  near-inertial  band.  Other  removal  mechanism  could 
be  the  wind  [Pollard  and  Millard,  1970],  and  in  the  next  section 
evidence  will  be  presented  that  shows  how  the  wind  by,  shifting  its 
direction,  may  remove  part  of  the  kinetic  energy  that  had  been  accumu¬ 
lated  earlier. 

Pollard  [1980]  shows  that  the  amplitude  of  inertial  oscillations  can 
fall  to  zero  within  a  few  tens  of  kilometers,  and  that  the  decay  of  energy 
is  due  to  a  slow  downward  dispersion  associated  with  horizontal  spreading, 
suplemented  by  the  effect  of  changing  winds.  In  the  present  work  there 
is  particular  interest  in  the  downward  dispersion  and  the  effect  of 
changing  wind  direction. 
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II.  THE  MILE  DATA 


A.  DESCRIPTION 

From  19  August  to  6  September  1977,  two  surface  moorings  were  de¬ 
ployed  in  the  vecinity  of  OWS  PAPA.  The  MILE-1  mooring,  which  provided 
the  data  used  in  this  study  was  located  at  49°  37'N  and  135°  6'W. 

This  mooring  consisted  of  an  8m  surface  toroid  anchored  with  0.9  scope 
in  4360m  of  water.  Shackled  in  the  mooring  line  were  19  VACM's  of 
which  the  upper  nine  were  used  in  this  study:  those  at  mean  depths  of 
5,  8,  11,  14,  20,  23,  26,  29  and  32m. 

The  reason  for  choosing  only  those  current  meters  was  that  the 
mixed  layer  was  at  all  times  more  shallow  than  about  30m,  so  the  measure¬ 
ments  below  32  meter  were  not  considered  to  be  of  much  significance  in 
testing  the  hypothesis. 

All  recorded  data  were  pre-processed  and  provided  on  magnetic  tape 
in  intervals  of  15  minutes  for  a  19-day  period,  with  the  following 
exceptions: 

At  the  5m  level,  there  are  only  5  days  of  current  record. 

At  the  8m  levle,  there  are  only  16  days  of  current  record. 

At  the  14m  level  there  are  only  3  days  of  current  record. 

Temperature  records  at  all  levels  are  complete. 

Wind  data  were  recorded  hourly  aboard  the  R/V  OCEANOGRAPHER. 

Figures  3  and  4  show  those  wind  values.  Two  strong  wind  events  can  be 
seen  one  strong  storm  in  22-24  August  (starting  time  100  hours),  and  a 
weaker  storm  in  the  31  August-1  September  period  (starting  time  310 
hours).  The  six-hourly  sea-surface  pressure  maps  reveal  that  both  storms 


Figure  4.  Win 


were  low  pressure  disturbances  whicli  developed  rapidly  in  the  Gulf  of 
Alaska  and  moved  toward  land  in  an  East-Northeast  direction.  The  center 
of  the  low  on  31  August  passed  directly  over  station  P,  causing  calm 
conditions  for  a  few  hours.  That  is  shown  by  the  double  pulse  in  figure  3. 

B.  ANALYSIS 

Figure  5  shows  hourly  temperatures  at  the  studied  levels  and  illus¬ 
trates  the  variability  of  the  temperature  structure  of  the  layer  above  the 
seasonal  thermocline.  Homogeneous  layers  are  indicated  in  this  figure  when 
temperature  traces  from  two  or  more  depths  coincide.  During  stratified 
conditions,  the  vertical  or  horizontal  movement  of  the  water  column  past 
the  sensors  causes  apparent  temperature  variability.  Most  of  this  vari¬ 
ability  appears  to  be  due  to  tidal-period  internal  waves.  Above  the 
30  meters  level,  homogeneous  and  stratified  conditions  alternate  in  re¬ 
sponse  to  atmospheric  conditions,  as  can  be  seen  in  the  power  spectra  in 
figure  5a. 

With  the  rapid  increase  of  the  wind  early  on  22  August,  at  about  95 
hours  on  the  time  axis,  the  temperature  at  5m  dropped  and  within  four 
hours  the  upper  23m  of  the  water  column  became  vertically  homogeneous 
due  to  strong  mixing.  During  this  strong  wind  event,  the  mixed  layer 
deepened  to  nearly  35m.  After  the  storm  died,  by  24  August  stratifi¬ 
cation  was  reestablished  between  17m  and  the  base  of  the  earlier  30m  mixed 
layer.  During  the  second  storm  the  weaker  wind  event  of  27-28  August 
mixed  the  upper  layers  to  about  22m  (figure  6). 

Figures  7,  8  and  9  are  the  power  spectra  for  the  u  and  v  com¬ 
ponents  of  the  observed  current  at  three  different  levels,  11,  23  and 


32m,  for  three  different  time  windows.  In  figure  7  the  time  window  is 
from  0  to  260  hours,  only  the  effects  of  the  first  storm  are  present. 

In  figure  9,  the  time  window  is  from  200  to  460  hours,  and  only  the 
second  storm  is  present.  For  the  spectrum  of  figure  8,  the  time  window 
is  from  100  to  360  hours,  and  the  time  period  includes  both  storms. 

Comparison  between  different  depth  levels  in  figures  7  and  9  indi¬ 
cates  different  amounts  of  kinetic  energy  for  each  depth  at  near-inertial 
frequencies.  The  first  storm,  as  can  be  seen  in  figure  4,  was  charac¬ 
terized  by  a  pronounced  shift  from  steady  easterly  winds  to  strong 
westerly  winds.  That  strong  change  after  a  period  of  at  least  100  hours 
of  steady  winds  could  be  the  means  by  which  part  of  the  energy  transfer 
between  the  ocean  and  the  atmosphere  could  have  eliminated  the  existing 
inertial  current  before  generation  of  a  new  inertial  current. 

Since  the  upper  ocean  is  assumed  to  be  one-dimensional  there  must 
be  some  unsteady  mechanism  by  which  the  vertical  transfer  of  momentum 
is  modulated  in  time,  causing  the  observed  phase  and  magnitude  differences 
between  the  32  meter  level  current  and  the  surface. 

The  response  of  the  current  in  the  upper  layers  to  the  second  storm 
was  quite  different  from  the  response  to  the  first  storm.  In  this  latter 
case  the  wind  shifted  from  having  a  northerly  component  to  a  strong 
southerly  one  and,  in  a  few  hours,  to  a  northerly  component  again 
(figure  4).  The  upper  level  currents  initially  appear  to  be  attenuated  by 
the  southern  component  winds  but  then  are  reinforced  when  the  wind  di¬ 
rection  shifts  toward  the  south.  Also  these  winds  have  a  smaller  inten¬ 
sity  and  duration,  and  as  it  is  shown  by  figure  6  the  mixed  layer  did  not 
deepen  as  much  as  in  the  first  storm. 
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Figure  7 


Power  spectrum  of  observed  current. 

Time  window  0  to  260  hours.  2 
(Abscissa:  c.p.h. ;  ordinate:  (cm/s)‘^) 
1  :  inertial  frequency. 
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Figure  9.  Power  spectrum  of  observed  current. 

Time  window  :  200  to  460  hours.  2 
(Abscissa:  c.p.h.;  ordinate:  (cm/s)  ) 
1  :  inertial  frequency. 
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Because  the  mixed  layer  had  undergone  a  shallowing  process  as  dis¬ 
cussed  earlier,  the  transfer  of  energy  to  the  32m  level  was  hindered  by 
stratification. 

Figures  10  and  11  are  the  power  spectra  for  the  u  and  v  compon¬ 
ents  of  the  current  for  the  total  period  of  observation.  The  spectra 
for  32m  level  exhibit  a  double-peaked  feature  at  the  near-inertial  period. 
As  Pollard  [1970]  suggests,  at  that  depth  it  is  unlikely  that  the  ampli¬ 
tude  of  inertial  oscillations  can  be  generated  by  the  action  of  a  single 
storm. 

The  time  series  of  observed  currents  are  shown  in  figures  12  to  19. 

In  those  figures,  the  plotted  running  mean,  an  average  over  one  inertial 
period,  shows  the  non-inertial  current  at  each  level.  It  should  be 
recognized  here,  that  an  undetermined  part  of  this  non-inertial  current 
may  be  quasigeostrophic  and  not  directly  related  to  the  local  wind 
conditions. 
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Figure  10.  Power  spectrum  of  observed  current. 
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(Abscissa:  c.p.h.;  ordinate:  (cm/s)  ) 
1  :  inertial  frequency. 
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Figure  15.  Observed  current:  0  to  260  hours  (v-component) . 

Running  mean  (*),  average  over  1  inert,  period. 
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Figure  17.  Observed  current:  100  to  360  hrs.  (v-component) 
Running  mean  (*),  average  over  1  inert,  period. 


III.  DESCRIPTION  OF  THE  MODEL 


A.  THEORETICAL  EQUATIONS 

In  this  study,  one-dimensional  conservation  of  momentum  in  the  mixed 
layer  is  assumed.  This  conservation  of  momentum  and  the  condition  of 
incompressibility  are  reflected  by  the  Navier-Stokes  equations  of  motion, 


Dv 

Dt 


+  7^ 


X  V 


-  4-  vp  +  S  + 


vV^v 


setting:  u  =  u^  +  u  +  u' 

V  =  Vg  +  V  +  v' 

where  u  and  u'  are  the  horizontal  mean  and  fluctuating  components  of 
x-component  of  non-geostrophic  velocity,  and  u^  is  its  geostrophic 
component;  similarly  v,  v'  and  v^  are  the  corresponding  y-components. 
With  the  Boussineq  approximation  and  Reynolds  averaging. 


9u'w'  9u' v'  u' u' 

9z  ’  9y  "  8x 


+  V 
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3v' v' 

'  3u' 

3y 

3x 

5  V 

+ 

3y^ 

Assuming 

(a)  Reynolds  stresses  are  much  greater  than  viscous  stresses. 

(b)  Vertical  velocity  field  vanishes,  w=0. 

(c)  Geostrophic  flow  is  separable  from  wind  driven  flow, 
and  neglecting  the  horizontal  density  gradient 


fu  =  - 

g  pQ  sy 


fv  = 

g  pQ  - 

Then  the  equations  for  wind-driven  flow  become 


=  _  fQ  _  v'w' 
3t  ^  3Z 


B.  NUMERICAL  METHOD  OF  SOLUTION 

To  solve  equations  [1]  and  [2]  a  vertical  grid  from  the  surface  to 
35m,  was  used.  The  values  of  u(z,t)  and  v(z,t)  were  then  solved 
numerically  as  an  initial-value  problem. 
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Taking  the  equations  for  wind-driven  flow  and  applying  an  eddy 
viscosity  closure. 


-u'w' 


-v'w' 


=  k 


3v 
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and  assuming  k  constant,  gives: 


3u  _ 
3t  ■ 

,2- 

fv  +  k  ^ 

3Z 

[3] 

3v 

3t  ■ 

-  fu  +  k  ^ 

3z 

[4] 

The  boundary  conditions  are: 
a.  At  the  surface: 


3u(0)  _ 
3z 


P 


-v'w^  = 


3v(0) 

3Z 


where  t  (t)  and  T„(t)  are  wind  stress  in  x  and  y 

A  y  '' 

directions  respectively  at  time  t  and  p  is  the  density 
of  sea  water. 

b.  At  z  =  -D  <  -h,  a  "slip"  condition  is  prescribed. 


u'w'(-D)  =  0 
v^(-D)  =  0 
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With  the  eddy  viscosity  closure. 


I  -  i  i-o)  -  » 

The  numerical  model  used  to  simulate  the  current  was  based  on  a 
"leapfrog"  scheme,  with  m  indicating  space  (depth)  and  n  time,  then 
u(z,t)  =  u(mAz,nAt)  =  u(m,n).  In  this  manner  the  component  of  equations 
[3]  and  [4]  can  be  written: 

au  u(m,n+l )  -  u(m,n-l ) 
at  "  2At 

av  _  v(m,n+l )  -  v(m,n-l ) 
at  ■  2At 

fv  =  fv(m,n) 
fu  =  fu(m,n) 

k  (u(m-l  ,n-l )  -  2u(m,n-l )  +  u{m+l  ,n-l ) ) 

az^  (Az)^ 

2 

k  (v(m-l  ,n-l )  -  2v(m,n-l )  +  v(m+l  ,n-l ) ) 

az  (Az) 

Solving  for  the  (n+1)  values, 

u(m,n+l )=Ru(m-l ,n-l )+(l-2R)u(m,n-l )+Ru{m+l ,n-l )+2fAtv(m,n) 
v(m,n+l )=Rv(m-l ,n-l )+(l-2R)v(m,n-l )+Rv(m+l ,n-l )-2fAtu(m,n) 
where:  R  =  2kAt/(Az)^. 

At  the  surface  m  =  M  and  application  of  the  boundary  conditions 
gives: 
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u(M,n+1)  =  2Ru(K-l.n-l)  +  (l-2R)u(M,n-l ) 

+  2fAtv(M,n)  +  2RAz  Tx(n)/(kp) 

v(M,n+l)  =  2Rv{M-1,n-1)  +  (1-2R)v(M,n-1 ) 

-  2fAtu(M,n)  +  2RA2Ty(n)/(kp) 

The  numerical  model  is  stable  when  the  value  of  R  is  equal  to  or 

less  than  0.5.  Assuming  an  eddy  viscosity  coefficient  k  of  23  cm/s 

and  an  increment  of  depth  of  Az  =  Im,  and  a  time  step  At  =  100 

seconds  gives  an  R  equal  to  0.46. 

The  initial  conditions  were  established  as  those  observed  values  of 

the  current  speed  for  each  component,  u  and  v,  at  time  0500  GKT 

August  19th,  at  the  first  recorded  peak  of  the  v-component  of  the  current. 

The  surface  stresses  tx  and  xy  were  computed  every  hour  as  a 

function  of  the  hourly  wind  speed  data  for  each  component  w  and  w  , 

X  y 

using 

2  2 

\  '  "x‘"x  *  "y> 

2  2  V2 

Ty  =  t  W^)  CpO, 

_3 

with  a  value  for  the  drag  coefficient  Cp  of  1.3  x  10  . 

Values  for  the  u  and  v  components  were  calculated  for  each 
level,  from  the  surface  to  35m  at  intervals  of  Im,  at  time  steps  of 
100  seconds. 

Every  hour  a  new  input  of  wind  stress  was  applied,  setting  as 
initial  conditions  for  the  next  run  an  integrated  value  for  depth  in  the 
mixed  layer. 
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Mixed  layer  depth  as  a  function  of  time  was  deduced  from  the  Mile-1 
temperature  time  series  and  checked  with  Plessey  CTD  data  taken  from  the 
R/V  OCEANOGRAPHER. 

The  model  was  initially  run  without  any  damping  term,  and  the 
results  v/ere  acceptable  until  the  onset  of  the  second  storm.  At  that 
time  a  clear  shift  in  phase  between  the  recorded  and  modeled  currents 
was  observed.  In  an  attempt  to  avoid  a  steadily  worsening  phase  problem 
a  linear  damping  term  was  introduced.  According  to  Pollard  and  Millard 
[1970],  this  decay  factor  models  the  dispersion  effect.  Then  the 
equations  for  wind  driven  flow  become; 

I  .  rt  .  I-  (UV,  .  ,5 

.  .rt  .  |_  (7n  -  cv 

-5  -1 

A  damping  coefficient  of  c  *  5.03  x  10  sec  was  chosen  from 
Pollard  [1970].  Thsi  gives  an  e-folding  decay  time  of  2.3  days. 
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IV.  RESULTS 


Figure  20  through  22  show  the  modeled  currents  compared  with  the 
recorded  ones  at  three  differents  levels. 

As  was  discussed  earlier,  the  dominant  frequencies  of  the  model 
and  the  ocean  are  slightly  different,  see  figures  23  and  24.  There¬ 
fore  when  the  wind  decreases  and  the  current  is  almost  totally  inertial 
under  small  forcing  the  modeled  and  observed  oscillations  will  drift 
out  of  phase.  This  is  most  apparent  in  the  period  before  the  second 
storm,  when  the  wind  speed  dropped  to  a  minimum. 

The  model  does  not  have  the  capability  to  reproduce  the  observed 
current  driven  by  the  sharp  change  in  directions  of  the  wind  during  the 
second  storm.  The  modeled  and  observed  currents  become  out  of  phase 
at  about  time  330  hours.  After  the  wind  again  becomes  more  steady  in 
direction,  the  modeled  and  observed  currents  again  agree  in  phase. 

Other  discrepancies  between  observed  and  modeled  current  could  be 
caused  either  by  observational  errors  or  by  failure  of  the  model  to 
properly  treat  the  physics  of  the  problem,  such  as  neglecting  advection. 

It  is  clear  that  the  dissipation  of  mean  kinetic  energy  is  too  much 
in  the  model  using  R  =  1/2,3  days.  Figure  22  shov/s  that  at  the  32m 
level  the  modeled  current  is  stronger  than  the  recorded  one. 

At  the  other  levels  the  opposite  happens.  The  model  needs  to  be 
tuned  to  get  a  better  amplitude  comparison  in  the  mixed  layer. 

Doing  this  tuning  by  only  comparison  with  observed  currents  is 
difficult  because  of  the  probable  errors  in  the  measurements.  As 
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Comparison  between  observed  and  modeled  currents 
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suggested  by  Hal  pern  et  al  [1981],  there  is  spurious  rotation  of  the 
rotor  of  the  current  meters  produced  by  motions  of  the  surface- 
following  buoy  mooring. 

Tuning  could  be  done  with  a  different  treatment  of  the  eddy  vis¬ 
cosity  as  well  as  alternate  parameterization  of  the  mean  kinetic  energy. 
The  first  approach  would  be  possible  when  turbulent  viscosity  of  the 
water  as  a  function  of  depth  and  stratification  is  better  known,  per¬ 
mitting  detailed  quantitative  calculations  of  the  vertical  stress 
profile. 

Two  other  factors  that  could  have  a  relative  important  effect  on  the 
currents  and  are  not  included  due  to  the  limited  scope  of  this  study  is 
the  process  of  inertial-gravitational  wave  propagation  and  the  resultant 
horizontal  dispersion  of  mean  kinetic  energy. 
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